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Abstract. Two particle correlations within a single jet produced in deeply inelastic scattering (DIS) off a 
large nucleus as well as in heavy-ion collisions are explored. This is performed within the framework of the 
medium modified dihadron fragmentation functions. The modification occurs due to gluon bremsstrahlung 
induced by multiple scattering. The modified fragmentation functions for dihadrons are found to follow 
closely that of single hadrons leading to a weak nuclear suppression of their ratios as measured by HERMES 
in DIS experiments. Meanwhile, a moderate medium enhancement of the near-side correlation of two high 
Pt hadrons is found in central heavy-ion collisions, partially due to trigger bias caused by the competition 
between parton energy loss and the initial Cronin effect. 

PACS. 12.38.Mh - ll.lO.Wx 


1 Introduction 

The modification of the properties of jets as they pass 
through dense matter has emerged as a new diagnostic 
tool in the study of the partonic structure of such an en¬ 
vironment Such modification goes beyond a mere sup¬ 
pression of the directed energy leading to a suppression of 
the multiplicity m and could in principle be extended 
to include the modification of many particle observables. 
Theoretically, such observables may be computed through 
a study of n-hadron fragmentation functions. Such func¬ 
tions may be defined as the n-hadron expectation values of 
partonic operators. These may be factorized from the hard 
collision and their evolution with energy scale and mod¬ 
ification in matter may be systematically studied in per¬ 
turbative Quantum Chromodynamics (pQCD) (20 ■ This 
is quite similar to the usual single inclusive hadron frag¬ 
mentation functions 

Experimentally, such objects may be estimated through 
the measurement of semi-inclusive distributions of hadrons 
emanating from such reactions. Along with the multiplic¬ 
ity distribution at a given momentum or momentum frac¬ 
tion, which corresponds to the single fragmentation func¬ 
tion, a variety of triggered distributions may also be mea¬ 
sured. In such experiments, the measurement of one or 
more associated hadrons is only undertaken in the event 
of the appearance of a trigger hadron with a certain char¬ 
acteristic momentum or momentum fraction. Such mea¬ 
surements have recently been performed in DIS off cold 
nuclear matter by the HERMES experiment at DESY |2]. 
Measurements in hot and dense nuclear matter have been 
performed at the Relativistic Heavy-Ion Collider (RHIC) 
by the STAR and PHENIX detectors |ni|7|. Due to issues 
related to background subtraction and requisite statistics, 


such measurements were restricted to two hadron corre¬ 
lations. While the two-hadron correlation is found to be 
slightly suppressed in DIS off a nucleus versus a nucleon 
target, it is moderately enhanced in central Au + Au colli¬ 
sions relative to that in p +p. This is in sharp contrast to 
the observed strong suppression of single inclusive spectra 
uni in both DIS and central A + A collisions. The theo¬ 
retical description of this behaviour forms the topic of this 
effort. 


2 Definitions 

Similar to single hadron fragmentation functions, dihadron 
fragmentation functions of a quark can be defined as the 
overlapping matrices of the quark fields and two-hadron 
final states. 
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and can be factorized from the hard processes (S), where 
p'^ = pf + P2 is the total momentum of the two hadrons 
with flavors hi and /12, 2 = 21 -|-22 is the corresponding to¬ 
tal momentum fraction and q± = pij_ —p2j_ is the relative 
transverse momentum. 

They also satisfy the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations as have been derived 
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in Ref. 0- One unique feature of the DGLAP equations 
for dihadron fragmentation functions is the contribution 
from independent fragmentation of two partons after the 
parton splitting, which involves the product of two single 
fragmentation functions. These DGLAP equations for di¬ 
hadron fragmentation functions have been solved numeri¬ 
cally m and the evolution of the dihadron fragmenta¬ 
tion functions agrees very well with results from JETSET 
CH Monte Carlo simulations of e+ -I- e hi + h 2 + X 
processes. Although both the single and dihadron frag¬ 
mentation functions evolve rapidly with their ratio 
has a very weak dependence. 

As in the case for the single fragmentation functions, 
dihadron fragmentation functions may not be predicted 
entirely within QCD, but need to be measured for all mo¬ 
mentum fractions at a given energy scale. Since there are 
no experimental data available for dihadron fragmentation 
functions, JETSET Monte Carlo results will be used as the 
initial condition for the vacuum dihadron fragmentation 
functions in this study. For single hadron fragmentation 
functions, the BKK parameterization m will be used; 
this also agrees well with JETSET results. 

Besides predicting the factorization and evolution of 
appropriately defined parton matrix elements, pQCD also 
predicts a universality of these matrix elements. The frag¬ 
mentation functions defined in e“'"e“ collisions are identical 
to that in Deep-Inelastic Scattering (DIS) and in hadron- 
hadron or nucleus-nucleus collisions. The hard scattering 
cross section ((Tq) is however different in each case. 


3 Modification in a cold medium 

Applying factorization to dihadron production in single 
jet events in DIS off a nucleus, e(Li) -|- A{p) e(L 2 ) + 
hi{pi) + h 2 {p 2 ) + X, one can obtain the dihadron semi- 
inclusive cross section as, 

ggpjg ^ a _^ avv 

(PL2dzidz2 27rs dzidz2 ’ 

in terms of the semi-inclusive tensor at leading twist, 
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q 
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In the above, (zi,Z 2 ) is the dihadron fragmenta¬ 
tion function, L^i, = (l/ 2 )Tr(/Li 7 ^ the factor 

represents the hard part of quark scattering with 
a virtual photon which carries a four-momentum q = 
[—Q^/2g“, 0j_] and fq{x, Q^) is the quark distribution 

in the nucleus which has a total momentum A[p'^, 0, 0^]. 
The hadron momentum fractions, zi = Pi /q~ and Z 2 = 
P 2 /q~, are defined with respect to the initial momentum 
q~ of the fragmenting quark. 

At next-to-leading twist, the dihadron semi-inclusive 
tensor receives contributions from multiple scattering of 


the struck quark off soft gluons inside the nucleus with in¬ 
duced gluon radiation. One can reorganize the total contri¬ 
bution (leading and next-to-leading twist) into a product 
of effective quark distribution in a nucleus, the hard part 
of photon-quark scattering and a modified dihadron 
fragmentation function 22 )- The calculation of 

the modified dihadron fragmentation function at the next- 
to-leading twist in a nucleus proceeds na similarly as that 
for the modified single hadron fragmentation functions [2] 
and yields. 
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In the above, xb = -Q‘^l2p+q , xl = l']_/2p+q y{l - y), 
l± is the transverse momentum of the radiated gluon, 
APq^qg and APq^gq are the modified splitting functions 
whose forms are identical to that in the modified single 
hadron fragmentation functions The switch {hi / 12 ) 
is only meant for the last term, which represents indepen¬ 
dent fragmentation of the quark and gluon after the in¬ 
duced bremsstrahlung. The corresponding modified split¬ 
ting function, 

.f, 1 + y^ CA2TrasT^g{xB,XL) 

l_y {ll + {kl))NJq^{xB,Q^)’ ^ ’ 

is similar to APq^qg but does not contain contributions 
from virtual corrections. In the above, Ca = 3, Nc = 3, 
and {k']_) is the average intrinsic parton transverse mo¬ 
mentum inside the nucleus. 

Note that both modified splitting functions depend 
on the quark-gluon correlation function in the nucleus 
T^, which also determines the modification of the sin¬ 
gle hadron fragmentation functions na For a Gaussian 
nuclear distribution, it can be estimated as |21 

T^g{xB,XL) = C{Q‘^)mNRAfq{xB){l - (6) 

where, xa = ^/ixinRa, totv is the nucleon mass and Ra = 
IA2A^A 

is the nuclear radius. The overall constant C oc 
x_lG(x_l) [a:_L = {k\)/2p'^q~y{l — y)] is the only parame¬ 
ter in the modified dihadron fragmentation function which 
might depend on the kinematics of the DIS process but 
is identical to the parameter in the modified single frag¬ 
mentation functions. In the phenomenological study of the 
single hadron fragmentation functions in DIS off nuclei, 
C = 0.006 GeV^ is determined within the kinematics of 
the HERMES experiment |H]. The predicted dependence 
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of nuclear modification of the single hadron fragmentation 
function on the momentum fraction z, initial quark energy 
u = q~ and the nuclear size Ra agrees very well with the 
HERMES experimental data |12j . 

With no additional parameters in Eq. 0 , one can pre¬ 
dict the nuclear modification of dihadron fragmentation 
functions within the same kinematics. Since dihadron frag¬ 
mentation functions are connected to single hadron frag¬ 
mentation functions via sum rules m, it is more illustra¬ 
tive to study the modification of the conditional distribu¬ 
tion for the second rank hadrons, 

R2{zi,Z2) = (7) 

where zi and Z 2 < zi are the momentum fractions of the 
triggered (leading) and associated (secondary) hadrons, 
respectively. Shown in Eig.QJis the predicted ratio of the 
associated hadron distribution in DIS off a nitrogen {A = 
14) target to that off a proton {A = 1), as compared to 
the HERMES experimental data 0. 

The modification of the dihadron fragmentation func¬ 
tion in a nucleus with atomic number A is estimated ex¬ 
perimentally by the HERMES experiment I3H1 by mea¬ 
suring the number of events with at least two hadrons with 
momentum fractions zi,Z 2 {N\{zi, Z 2 ))i and the number 
of events with at least one hadron with momentum frac¬ 
tion z {Na{z)). In an effort to eradicate systematic errors, 
a double ratio is presented (using similar measurements 
off deuterium, i.e., A = 2), 


”^2(^2) 


S., >0.5 


( 8 ) 


This double ratio for different Z 2 is plotted in Fig.^as the 
circular points.. 

In comparison with the experiment we present the as¬ 
sociated hadron distribution: the ratio of the modified di¬ 
hadron to the single hadron fragmentation function in a 
nucleus normalized by the same ratio in the vacuum. As 
in the experiment, we integrate over zi > 0.5 in both the 
single and dihadron fragmentation functions. The agree¬ 
ment between the prediction and the data is remarkable 
given that no free parameters are used. The suppression 
of the associated hadron distribution at large Z 2 due to 
multiple scattering and induced gluon bremsstrahlung in 
a nucleus is quite small compared to the suppression of the 
single fragmentation functions iHiini. Since 77-2 is the ratio 
of double and single hadron fragmentation functions, the 
effect of induced gluon radiation or quark energy loss is 
mainly borne by the single spectra of the leading hadrons. 
This is similar to the evolution of the dihadron fragmen¬ 
tation function with momentum scale in the vacuum |31 
n Oj . At small values of Z 2 , the modified dihadron fragmen¬ 
tation function rises above its vacuum counterpart more 
than the single fragmentation functions. This is due to the 
new contribution where each of the detected hadrons em¬ 
anates from the independent fragmentation of the quark 
and the radiated gluon. 



Fig. 1. Results of the medium modification of the associ¬ 
ated hadron distribution, in a cold nuclear medium versus its 
momentum fraction. The momentum fraction of the leading 
hadron zi is integrated over all allowed values above 0.5. The 
ratio of this quantity in Nitrogen (N) with that in deuterium 
(D) are the HERMES data points. 


4 Modification in a Hot medium 

In high-energy heavy-ion (or p + p and p A) collisions, 
jets are always produced in back-to-back pairs. Correla¬ 
tions of two high-pT hadrons in azimuthal angle generally 
have two Gaussian peaks EE Relative to the triggered 
hadron, away-side hadrons come from the fragmentation 
of the away-side jet and are related to single hadron frag¬ 
mentation functions. On the other hand, near-side hadrons 
come from the fragmentation of the same jet as the trig¬ 
gered hadron and therefore are related to dihadron frag¬ 
mentation functions. 

To extend the study of intrajet dihadron correlations 
to heavy-ion collisions, one can simply assume {k\) ~ /i^ 
(the Debye screening mass) and that C in Eq. 0 is pro¬ 
portional to the local gluon density of the produced dense 
matter. In addition, one also has to include the effect of 
thermal gluon absorption CSl such that the effective en¬ 
ergy dependence of the energy loss will be different. Such a 
procedure was applied to the study of the modification of 
the single hadron fragmentation functions and it success¬ 
fully described the quenching of single inclusive hadron 
spectra, their azimuthal anisotropy and suppression of 
away-side high pT hadron correlations m The change 
of the near-side correlation due to the modification of di¬ 
hadron fragmentation functions in heavy-ion collisions can 
be similarly calculated. 

The near-side correlation with background subtracted 
and integrated over the azimuthal angle m can be re¬ 
lated to the associated hadron distribution or the ratio 
of dihadron to the single hadron fragmentation functions 
as in Eq. O- However, one has to average over the ini¬ 
tial jet energy weighted with the corresponding production 
cross sections. For a given value of of the triggered 
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hadron, one can calculate the average initial jet energy 
{Et)- Because of trigger bias and parton energy loss, {Et) 
in heavy-ion collisions is generally larger than that in p+p 
collisions for a fixed In these proceedings, we simply 
use the calculated {Et) and its centrality dependence for 
a given from Ref. CHI which in turn determines the 
average value of zi = p)p^/{E t) and Z 2 = Pt^°^/{Et)- 
In the same calculation, one can also obtain the average 
parton energy loss as the difference of the average ini¬ 
tial parton energy with and without parton energy loss, 
{AE) = which is completely de¬ 

termined from the suppression of single inclusive hadron 
spectra and away-side correlation. Since the high bi¬ 
ases the jet production position towards the surface of 
the overlapped region, {AE) associated with a triggered 
hadron is always smaller than the average energy loss of 
both the away-side jet and jets in non-triggered events. 

Using {Az) = {AE)/{Et), we first determine the pa¬ 
rameter C and in turn calculate both the modified single 
and dihadron distributions. The ratio of such associated 
hadron distributions in Au + Au versus p + p collisions, 
referred to as Iaa [H]; is plotted as the solid line in Fig.0 
together with the STAR data [3, as a function of the 
number of participant nucleons. In central Au + Au col¬ 
lisions, triggering on a high pT hadron biases toward a 
larger initial jet energy and therefore smaller zi and Z 2 . 
This leads to an enhancement in Iaa due to the shape 
of the dihadron fragmentation functions P]. This is in 
contrast to the slight nuclear suppression of the associ¬ 
ated hadron distribution at large Z 2 in DIS off a nucleus 
(Fig. [TJ. The enhancement increases with Npart because 
of increased initial gluon density and system size of the 
produced dense matter which leads to an increased total 
energy loss. In the most peripheral collisions, the effect 
of smaller energy loss is countered by the Cronin effect 
due to initial state multiple scattering that biases toward 
smaller {Et) relative to p -|-p collisions. As a result, the 
associated hadron distribution is slightly suppressed. 


5 Conclusions 


In summary, we have studied the medium modihcation 
of dihadron fragmentation functions in both hot and nu¬ 
clear matter due to multiple parton scattering and induced 
gluon radiation. The modification is found to follow closely 
that of single hadron fragmentation functions so that the 
associated hadron distributions or the ratios of dihadron 
and single fragmentation functions are only slightly sup¬ 
pressed in DIS off nuclei but enhanced in central heavy-ion 
collisions due to trigger bias. With no extra parameters, 
our calculations agree very well with experimental data. 

This work was supported by the U.S. Department of 
Energy under Contract No. DE-AC03-76SF00098, NSEC 
under project Nos. 10475031 and 10135030 and NSERC 
of Canada. 



Fig. 2. Calculated medium modification of associated hadron 
distribntion from jet fragmentation in Au -|- Au collisions at 
Vs = 200 GeV. E]. 
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